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ABSTRACT 


The results of tests made with a two dimensional diffusor 
operating with transverse grooves cut into the diverging side walls 
were compared with the results of similar tests made using a smooth 
wall diffusor. The results of this comparison revealed that: 

1. The efficiency of the ribbed wall configuration exceeded that 
of the smooth wall diffusor for total opening angles of approximately 
28 to 33 degrees. 

2. The ribbed wall configuration had an efficiency approximately 
120 per cent of that of the smooth wall diffusor when both were oper- 
ated at a total opening angle of 32 degrees. 

3. The efficiencies of the ribbed wall diffusor operating at 
large opening angles varied as a function of the number of grooves in 
the side walls. 

4, The efficiency increased as the number of grooves increased, 
but an unexplained drop in efficiency occurred at openitig angles of 
30 and 32 degrees with an eight groove configuration. 

5« The maximum number of grooves tested was ten. There is no 


evidence to indicate that this is an optimum configuration. 
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CHAPTER I 
INTRODUCTION 


Much work has been done and many papers written concerning studies 
of diffusor performance. Most of these studies are limited to the 
study of diffusors operating with a total opening angle of less than 
15 degrees. Cocanower, Kline, and Johnston, Ref. ', have proposed a 
general theory for predicting diffusor performance for a variety of 
two and three dimensional diffusors. Their theory has been shown to 
be in good agreement with experimental tests. 

The studies of diffusor performance for diffusors operating at 
large opening angles have not progressed as rapidly as the small angle 
diffusor work because of the complications associated with the severe 
turbulence resulting when flow separation occurs at the diffusor walls. 
This separation of flow is responsible for the large losses and poor 
performance normally associated with wide angle diffusors. 

Two studies which deal with the problem of wide angle diffusors 
are the papers of Waitman, Reneau, and Kline, Ref. 2, and of Moore, 
Carl, and Kline, Ref. 3. Each of these studies deals with the condi- 
tions affecting flow separation in wide angle diffusors. 

The flow separation phenomena observed and discussed by Moore, 
Carl, and Kline in connection with water table tests were also observed 
in the air flow tests of this investigation. These phenomena are dis- 
cussed in Section 3.6 of Chapter III. 

Migay, Ref. 4, noted an improvement in the efficiencies of dif- 
fusors operating at wide opening angles when transverse grooves were 


cut into the diverging walls. His explanation for this improvement in 


efficiency was that the ''large separation'' normally occurring without 
the grooves is replaced, because of their effect, by a series of 
‘smaller separations.'' The resulting losses associated with these 
"smaller separations'' are apparently less than the loss associated with 
the "large separation.'' The efficiency is thereby improved. 

Migay cites the results of many tests which support his claim of 
an improved efficiency when operating with a grooved configuration. 
Unfortunately, many of the results cited are not adequately documented. 
The results of the two dimensional experimentation, for example, are 
presented as the ratio of ribbed wall diffusor efficiency to that of a 
smooth wall diffusor operating under the same conditions. No informa- 
tion 1s available concerning the specific conditions under which the 
tests were conducted, and the actual efficiencies of each of the dif- 
fusors are not presented in the paper. 

The over-all objective of this study was to determine experimen- 
tally the performance of a two dimensional ribbed wall diffusor. It 
was hoped that the optimum operating range of such a diffusor could be 
found and compared with a smooth walled diffusor operating under simi- 
lar conditions. 

Figure 1 illustrates the type of diffusor used for the experiment. 
The detailed insert in the figure shows what is meant by the term 
"ribbed wall.'! 

Because of time limitations, the scope of this study was narrowed 
to the study of one particular rib configuration. The geometry of 
this configuration is discussed in Section 2.4 of Chapter II. A limited 
investigation was also made to determine the effect that varying the 


number of grooves had upon the performance of the ribbed wall diffusor. 


Liars 


TWO DIMENSIONAL DIFFUSOR 


ITEM A 


DETAIL OF RIBBED WALL CONFIGURATION 


FIGURE | 
TWO DIMENSIONAL DIFFUSOR AND 


DETAIL OF A RIBBED WALL 
CONFIGURATION 
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Migay cites a limited amount of data related to a two dimensional 
ribbed wall diffusor. The data were obtained using a one sided dif- 
fusor with an asymmetric opening angle of 20 degrees. Figure 2 shows 
the type of diffusor used. Table I is a summary of the information 
presented in Ref. 4 related to two dimensional diffusors. 

The initial rib configuration utilized for this study was chosen 
to be approximately equal to that used in Ref. 1. In this study, how- 


ever, a symmetric diffusor was used. 
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TABLE I 


SUMMARY OF TWO DIMENSIONAL DIFFUSOR STUDIES FROM REF. 1 


Effect of Parameter A (B = 5mm 


A/B 0.0 Calor Oe250 


CO 
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CHAPTER II 
EXPERIMENTAL PROCEDURES 


2.1 Experimental set-up. The equipment used for this study was 
located in the ''New Addition'' of Building 230 at the Jet Propulsion 
Laboratories, Naval Postgraduate School, Monterey, California. Figure 
3 1S a schematic of the installation. The source of air for the tests 
was the Carrier Air Compressor installed in Building 230. 

2.2 Pressure readings. All pressure measurements in the flow 
Field were made using a Fused Quartz Precision Pressure Gage manuf ac- 
tured by Texas Instruments Inc. of Houston, Texas. Prior to, and twice 
during the time period of the tests, the pressure instrument was cali- 
brated using a water manometer. The calibration remained invariant 
throughout the duration of the test; and the value of 0.8 inches of 
water pressure per increment, or count, on the instrument was used as 
the conversion coefficient for all pressure measurements. The details 
of the conversion of the pressure and temperature readings are dis- 
cussed in Chapter III, Sections 3.1 and 3.2. 

Atmospheric pressures were read from a standard wall barometer. 
Two pressure readings were required at the flow nozzle. The static 


pressure upstream, P,,, and the drop in static pressure across the 


u 
nozzle, 4P , were measured and recorded for the purpose of determining 
the flow rate through the system. The computations involved in the 
determination of the flow rate are discussed in Chapter III, Section 

5 ISI 


Two pressures were measured at the inlet plane of the diffusor, 


the static pressure and the dynamic head, a » The static pressure 
4 
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was assumed to be uniform in the plane of the throat. Two static taps 
were located at the top of the diffusor inlet and two at the bottom. 
All four of these static pressure taps were connected to the pressure 
measuring instrument through a common lead. The total pressure was 
measured at three positions in the central region of the flow, free 

of the boundary layer effect. Figure 4 shows the location of the pres- 
sure orifices and the details of the total pressure probe construction. 
The probe was manufactured by Mr. M. 0 dell, machinist, Aeronautics 
Department, NPGS. 

Figure 5 is a schematic of the pressure instrument hook-up. By 
opening and closing various pairs of the six switches, it was possible 
to measure all of the required pressures. 

Section 3.4 of Chapter III illustrates how the conditions at the 
throat were calculated using these pressure measurements. 

2.3 Temperature measurements. Two total temperature readings 
were taken in the flow field. The first was taken just ahead of the 
Flow nozzle and was used in the flow rate calculation. The second 
total temperature measurement was taken just ahead of the diffusor 
inlet. This temperature was used in the determination of the flow 
conditions at the inlet. 

All flow field temperatures were measured using a thermocouple. 
Section 3.1 of Chapter III discusses the details of the temperature 
reading conversions. The ambient air temperatures were measured with 
ordinary mercury bulb thermometers. 

2.4 Measuring and recording the data. The rate of flow through 


the system was regulated by bleeding the compressor output to the 
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FIGURE 4& 


SCHEMATIC OF THE DIFFUSOR INLET PLANE ILLUSTRATING THE LOCATION 
OF THE TOTAL AND STATIC PRESSURE TAPS 
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atmosphere. The maximum obtainable flow rate was approximately 2.4 
Ib /sec. At this flow rate, the velocity in the center of the inlet 
flow channel was approximately 210 ft/sec. 

In the initial phase of the experiment tests were made at diffusor 
angles varying from /3 = 10° to /3 = 60°. At each angle, the flow rate 
was varied between the maximum of 2.4 and a minimum of 0.8 1b,/sec. 

Tests were made at each of the angles with smooth walls and with 


ribbed walls configured as indicated below: 


Ns: number of grooves = 10 

b: width of grooves = 0, 125% 
t: spacing between grooves = 0.060 in. 
a: depth of grooves = 0.700 Gn. 


This rib configuration was chosen as a close approximation to the rib 
sizes used by Migay in the two dimensional diffusor study cited earlier 
and summarized in Table I. 

The results of this phase of the experiment indicated that the 
ribbed walls performed best, relative to the smooth walls, at [> =g20 
and a velocity of approximately 155 ft/sec. Consequently, further 
tests were conducted at this operating point in an attempt to determine 
how the performance varied as a function of the number of grooves. In 
this phese of the experiment the number of grooves was increased by 
one for each successive run. The results of these tests are discussed 
in Chapter IV. 

The program used to perform the required calculations is discussed 
in detail in Appendix A. A sample data sheet, a list of the variable 
names used on the sample data sheet, and a listing of input variables 


are included with the discussion of the program. 


Ze 


CHAPTER III 
COMPUTATIONS AND ANALYSIS 


3.1 Temperature measurement conversions. The thermocouple 
readings, recorded in millivolts, were converted to degrees Fahrenheit 
by the following equation: 

PF = B24 35.9B8(mv) -. 435 (my) * 
(3-1) 

This equation is a least squares fit of the values for conversion 
from millivolts to degrees Fahrenheit as specified in the National 
Bureau of Standards Circular 561 for iron constantan thermocouples. 
Equation (3-1) is valid only if the reference junction is maintained 
at a temperature of 32 degrees Fahrenheit. 

3.2 Pressure measurement conversions. The barometric pressure, 
measured and recorded in inches of mercury, was converted to pounds per 
Square inch by using the following equation: 

0001634 ¢t 
Pt 9120 6 
I+ 0001818 t 
(3-2) 
The two coefficients multiplying t, the temperature in degrees centi- 
grade, are thermal expansion coefficients for the mercury in the baro- 
meter and the brass upon which the scale of the barometer is engraved. 
The factor .4912 is the conversion factor for converting pressure 
measured in units of inches of mercury to units of lbs. per sq. in. 
Equation (3-1), the conversion factor, and the correction coefficients 
were obtained from Ref. 5, pages 1-157. 
3-3 Flow rate measurements. References 6 and 7 provide general 


methods, procedures, and calculations concerning flow rate measurements 


23 


using flow metering orifices and nozzles. The flow metering nozzle 
used to determine the flow rate for these tests was constructed in 
accordance with the specifications of Ref. 8. A nozzle was selected, 
rather than a square-edged orifice, because the calculations were felt 
to be more easily adaptable to computerized methods. 

From Ref. 8, page 5, Eq. (la), slightly modified, the flow rate 


through the flow metering nozzle in 1b,/sec is: 
wes Ae Faas (2 g. o., (hu- Fa) 


The quantities in Eq. (3-3) are identified below: 


(3-3) 

v/s flow rate (1b,,/sec) 
O% discharge coefficient 
é density coefficient 
A\ area of nozzle opening (in.“) 
Co area multiplier 

£. Reynolds number correction 

a. conversion factor (32.17 ft. lb,,/sec* Ibe) 

(R,- Pa) pressure drop across nozzle (psi) 

by upstream ''density'' (non-standard units) 

By assuming the flow to be a perfect gas, d.,08" be expressed as 


a function of the temperature, Tu. 


Pu 
NT 
a (3-4) 
Ra gas constant for air (53.35 feb 


Tu temperature upstream (°R) 
The "equivalent flow rate,'! denoted as W » and defined by Eq. 


(3-5), provided a computational means of determining the flow rate. 
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Pa QV 
25) 


& 
In Eq. (3-5),W has the units of ins’. 


After substituting for We and Ce, Eq. (3-5) becomes: 


| “A Ze: 


The table below lists the values of the coefficients in Eq. (3-6) 


which are applicable to the particular nozzle installed for these tests. 


TABLE II 


SUMMARY OF FLOW RATE CALCULATION COEFFICIENTS 


Coefficient Symbol Coefficient Value Source of Value 





ig xs Fash » then the coefficient E& can be expressed by the 


equation, 





vue ee a9 re Lun 
| x \- B05 (j- xy" 


(3-7) 

If t, is the temperature of the flow upstream of the nozzle, ex- 
pressed in degrees Fahrenheit, then the coefficient 1S given by the 
equation 

bar ee) 


f, = 1+ 001 (S852= 
lOO (3-8) 


a5 





Equation (3-6) may be greatly simplified by combining all of the 


known constant coefficients into one constant, c. 


C= x fF, Ay |2 =\\.30 


Thus, Eq. (3-6) becomes: 
= 11.30 € Fy] Cos) 
(3-9) 


Finally, the mass flow rate can be calculated by: 


ve RA 


= {Ta | R (3-10) 


3.4 Flow conditions at the diffusor inlet. Velocity and total 
pressure measurements across the inlet of the diffusor indicated that 
the static pressure was uniform throughout the inlet plane. It was 
assumed that the total pressure at the inlet was constant in the region 
outside of the boundary layer. Based upon this assumption, phe Meer 


number in the central region of the flow at the inlet was determined 


by solving the following isentropic relationship for My: 


Beef te cGaeante nm") +I 
P 2 , (3-11) 


The measured values of 7 and Re were used to determine My ° 


set 1 0G he | ae 


The static temperature at the inlet was determined by using the 


previously determined Mach number and the following equation: 


(10 em) 


(3413) 
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Solving Eq. (3-13) for TV, yields the required relationship: 


a MM (3-14) 
Once il » the static temperature, was determined, it was possible 
to evaluate a, the speed of sound, and subsequently the micstream velo- 
city at the inlet. 
— oa MO, - mM, Ro GoM 
(3-15) 
The velocity in the central region of the flow is indicated by V,.,- 
By assuming the flow to be a perfect gas, enough information was 
available to calculate the density of the flow at the inlet. It was 
assumed that the density remained constant across the inlet plane. 
Qe - Pr x 144 
Ra Be Ui (3-16) 
Equation (3-16) expresses the density in units of slugs/ft3. 
The Reynolds number based upon the equivalent diameter of the 
inlet was determined by expressing the kinematic viscosity of air as 
a function of temperature. The expression represents a fit of the 
data tabulated on page 8 of Ref. 9. The following expression is con 
sidered valid only for the temperature range of 50°F. to 104°F 
D,= (0.6 t+ 120)* \O7° 
(3-17) 
The units of 2) are ft*/sec. 
3.5 Integrated flow rate at the diffusor inlet. All of the cal- 
culations related to the conditions at the inlet and the performance 


of the diffusor were dependent upon the total pressure, static pressure, 


and total temperature measurements at the diffusor inlet. It was felt 
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that some method of evaluating the reliability of these measurements 
should be incorporated into the computational procedures. The method 
decided upon was to compare the flow rate at the inlet, computed from 
the velocity and density as determined from these measurements, against 
the flow rate measured at the flow nozzle. Since the two flow rates 
were determined from completely independent measurements, it was felt 
that any erroneous readings, or errors in recording the data, would 
become apparent when the two flow rates were compared. 

Prior to performing the calculations to determine the flow rate 
through the diffusor, tests were made with a hot wire anemometer to 
determine the velocity profiles at the inlet. Figure 6 is the velocity 
ratio, V/Vmax, plotted as a function of the distance, h, from the sur- 
face of the wall. Vp,, is the velocity in the central region of the 
Flow. Figure 7 represents the same velocity ratio as measured from the 
side wall. The profiles represented in Figs. 6 and 7 were measured at 
the mid-span position of the noe and side walls respectively. 

In each instance, it was possible to find analytic expressions 
which closely approximated the velocity profiles. The data plotted 
were taken at several velocities representative of all anticipated 
operating conditions. Although it was recognized that the velocity 
profiles vary as a function of velocity, tests indicated that the vari- 
ation was minor and under the anticipated velocity operating range 
would not exceed the variations caused solely by instrumentation error. 
The total change in the flow rate determination over the extremes in 
the velocity profile predictions was found to be a fraction. of 1 per 


cent of the total flow rate. 
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The expression used to approximate the velocity along the top wall 


- 7 160 ( 1 Vere 


h= 120 € 


(3-18) 
When solved for V as a function of h, the distance from the wall, 


Eq. (3-18) becomes: 
’ Sole 1.2O0\A 


A similar expression was obtained for the shape of the velocity 


(3-19) 


profile at the side wall. 


——,.7 97130 Cl- Wirax 


(3-20) 
Solving for V as a function of h yields: 
: le 0.60 \4 
V = Vian {| 30 (dn aot | (3-21) 


Figure 8 is a cross-sectional view of the inlet plane. The region 
labeled ''C'' is the area outside of the effect of the boundary layer. 

It was assumed that the velocity in this region, previously referred to 
as Vee was constant throughout. Measurements indicated that this 
assumption was reasonable. 

A turbulence survey was made in the central region by feeding the 
hot wire anemometer output to an oscilloscope. The results of this 
study indicated that the 1 maximum velocity perturbation was less than 
Wiperecentwof thesmean velocity, 

The regions labeled '"'F'' and ''D'' in Fig. 8 represent those regions 
within the boundary layer along the top, bottom, and side walls 


respectively. 


cw 


JFINI 3JHL HONOYHL J1LVY MOT1S JHL 3O NOTIVNIWY3130 JHL NI G3SN SV 
SNOT9ISY MOTI ONY SYFAVT AYVONNOG JO NOILVIO1 JHL SNIMOHS YOSNSSIO 4O 3NV1d LIINI 


8 Jynoi4 


19° 0 nQ°r 19°O 





10° 4 


10°9 


Sy 


Although it was not possible to accurately measure the velocity 
profiles on the bottom and opposite side wall, those measurements taken 
indicated that the variations from wall to wall and top to bottom were 
slight. Thus, it was assumed that the velocity profiles on opposite 
walls were identical. 

The corners, indicated as ''E'' in Fig. 8, were assumed to have a 
velocity distribution equal to the product of the W/V ratios from 
each wall as a function of the respective distances. Figure 9 shows 
how this was done. The velocity assumed to exist at the point p, 
indicatca if, Fige 9, was assumed to be the product of the velocity 
ratio ae given by Eq. (3-19) for the distance h,, from the top 


wall and the velocity ratio V/V given by Eq. (3-20), for the 


max? 
specified distance, hy, from the side wall. 
The diagonal line indicated as "corner survey path line'' in 
° 
Fig. 9 is the line upon which a survey of velocities was made. Figure 
10 shows the comparison between what the above assumption predicted for 
the velocities along the line and the actual measured values of V/V,,,y- 


The analytic expression for the velocity at any point in the corner 


region is given as a function of the two distances h, and ho. 


V2 Var lee bn EY ELI BoC kn SEF a 


The purpose of this investigation into the velocity profile at the 
inlet was to obtain a measure of the reliability, or validity, of the 
pressure readings. This was accomplished by using the measured values 
of the velocity and density together with the estimates of the velocity 


profile to integrate across the inlet and determine the flow rate. 
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CORNER SURVEY PATH LINE 


FIGURE 9 


THE CORNER REGION (E) OF THE INLET PLANE 
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The calculation for the mass flow rate is expressed by the fol low- 


ing equation when the flow is considered to be incompressible: 


S| 
v8 Sgeb Vv 4h ans 


Table III is a breakdown of the mass flow rate calculations for 


each of the regions specified in Fig. 9. 


TABLE III 


SUMMARY OF CONTINUITY CALCULATIONS 


Region Mass Flow CT SEC ) | Eqs 


G3. Vmax A 
0.6 
Eq. (3-19) 0.022 04. § V. dh, 


127 oy 


Eq. (3-22)! 0,028 eg. ( \ V, aahydh, 
Eq. (3-21)| 0,067 eg. * Vz ane 





The flow rate determined by this integration was compared to the 
Flow rate measured at the flow nozzle as discussed in Section 3.2. 

In general, the agreement was within 1 per cent. 

3.6 Diffusor performance. The conditions which existed at the 
diffusor inlet were found to be generally steady, and the measurements 
were felt to be reliable indicators of the flow at that plane. It was 
found, however, that conditions at the exit plane were, under most 
operating conditions, extremely turbulent. For large opening angles, 
of particular interest in this study, the flow actually reversed direc- 
tion along the diffusor walls after separation had occurred. Moreover, 


at opening angles greater than 26 degrees, the flow was unsteady in that 
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it would separate first on one side of the diffusor and then on the 
other. It was noted that under certain conditions the flow could be 
directed along one wall of the diffusor while complete separation ex- 
isted along the opposite wall. The flow could then be shifted to the 
opposite wall causing the separation to shift sides. It was also noted 
that the diffusor could be made to separate at both sides simultaneously 
causing a ''free jet'' flow straight out of the inlet throat. 

Because of these non-uniform conditions which existed at the exit, 
it was not possible to measure a single meaningful exit velocity. The 
only flow condition which could be accurately measured at the exit was 
the static pressure, Py. Because the diffusor exhausted to atmospheric 
conditions at subsonic velocities, the exit static pressure was taken 
to be equal to atmospheric pressure. 

Not being able to measure an exit velocity posed a serious problem 
in attempting to measure diffusor performance. Most definitions of 
diffusor efficiency, for example, require a knowledge of the exit velo- 


eye, VO. Wigay, Ref. |, defines wdiitfuser erficiency as: 


YY - Pe a P 
> HR(V%V,") (3-28) 
He makes no attempt to explain the meaning of Vo or how it was obtained. 
A more convenient and perhaps realistic criteria for diffusor per- 


formance are the coefficient of pressure recovery, defined as 


o* al P2 -% = 
* ARV, (3-29) 


This measure of diffusor performance has the distinct advantage 
that all of the quantities can be easily determined. Vms,, the velo- 
city of the flow in ''Region C,'' Fig. 8, was chosen as the reference 
Velocity Vj, to be used in Eq. (3-29). 


By) 


In order to compare the results of this test with those published 
in Ref. 1, some scheme had to be decided upon to produce a quantity V9 
for Eq. (3-28). One method used was simply to assume that the central 
flow at the inlet, with velocity V,.,, was expanded completely and in- 
compressibly from area Aj to area Aj. Vo is an ideal expansion of the 
midstream velocity at the inlet. 

A second means of evaluating the efficiency was to determine the 
ratio of the measured to the ideal coefficient of pressure recovery. 

The ideal coefficient of pressure recovery can be expressed as a 


function of the inlet and exit areas, 


Cpr, 7 -(M) 


Equation (3-30) was derived by substituting Bernouilli's equation 


(3-30) 


and the equation of continuity for incompressible flow into Eq. (3-29). 
3.7 Measurement reliability indicator, In order to obtain an- 
other indication of the reliability of the pressure and temperature 
measurements at the diffusor inlet, a comparison was made of the mea- 
sured flow rate based upon the flow metering nozzle and the flow rate 
calculated by Pradr ating the velocity over the area of the diffusor 
inlet as described in Section 3.4. The quantity referred to as the 


"data reliability indicator'’ was defined by the equation: 


ERR =| Wnioa 
Winter’ (223i) 
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CHAPTER IV 
RESULTS AND RECOMMENDATIONS 


4.1 Results. This project was conducted in three phases. The 
preliminary phase consisted of planning, building and evaluating the 
test equipment. The second phase consisted of the development of 
procedures and data reduction calculations. The final phase encom- 
passed the testing required to determine the most desirable operating 
range for conducting further tests and an investigation to determine 
the effect of the number of grooves upon the ribbed wall diffusor per- 
formance. The results of each of these phases of the project are dis- 
cussed individually in the following sections. 

4.1.1 Preliminary phase. The comparison of the flow rate meas- 
ured at the nozzle and the flow rate computed at the diffusor inlet 
provided a cross check of the pressure readings and the computed 
values of the Mach number and velocity. The data reliability indica- 
tor defined by Eq. (3-31) reflected the comparison of these two flow 
rates. The average value of this indicator for the last 50 data 
points recorded during the project was 0.0044, This indicates that 
the agreement was, on the average, within 0.5 per cent. This is com- 
parable with the accuracy of the flow nozzle which is * 0.5 per cent. 

Section 3.6 of chapter III discusses the two methods used to 
determine the diffusor efficiency. [It can be shown that the ratio 
PR /cog, is identical to the efficiency defined by Eq. (3-28). Both 
methods of determining the efficiency were included in the calculations 
dealing with diffusor performance since Eq. (3-28) is based upon the 


calculated value of the density, @ , whereas the ratio of pressure 


BI 

















recovery coefficients does not depend upon ‘4 » the comparison of 
these two efficiencies provided, therefore, an indication of the 
reliability of the measured value of the temperature and pressures at 
the inlet. 

The third indication of the reliability of the test equipment 
and procedures was the repeatability of the tests. Table IV shows 
the results of tests run on separate days under similar operating 
conditions. It was felt that the repeatability indicated in Table IV 
was reasonably good and that a realistic performance determination 
could be made on the basis of these tests. 

The results of the tests seemed to indicate that the test equip- 
ment and procedures could be relied upon to an extent consistent 
with the intended use of the calculated results. It was felt that 
sufficient safeguards had been incorporated into the computational 
system to provide adequate warning indications in the event that 
erroneous data were being processed, 

4.1.2 Determination of the most desirable operating range. To 
accomplish the objective of this phase of the study, ribbed wall and 
smooth wall tests were run over a wide range of operating conditions. 

Table V summarizes the results of this phase of the tests. Al- 
though runs were made from f3 10° te B = 60° for velocities 
ranging from 65 ft/sec to 210 ft/sec, only those values of (3 which 
were of particular interest to this study are presented in Table IV. 

Figure 11 illustrates the operating range of particular interest 
for ribbed wall investigations. Figure 12 shows how the ribbed wall 


efficiency compares to that of the smooth wall diffusor in this region. 
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TABLE IV 


RESULTS OF REPETITIVE TESTS UNDER SIMILAR 
OPERATING CONDITIONS 


Ribbed Wall Smooth Wall 
Tne: [Mme | mom | Oy 
220 













B 
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TABLE V 


COMPARISON OF SMOOTH WALL AND RIBBED WALL 
DIFFUSOR EFFICIENCIES 


S3200met ---—— | = 
22° Sipe 78.29 0.9640 
he ee hl ees 69.79 | ------ 
26° 65.30 elo 0.9456 
28° 58.58 59.63 1.0179 
30° BS joe | 58 uk 1.1015 
32° 47.84 56.95 1.1906 
34° 43.51 34.55 0.7941 
36° 39.30 29.82 0.7588 
380 3 (enc —— ee | ae 
40° SI otele 27.69 0.8746 


*All data recorded at a measured flow rate of 
approximately 1.78 1b,/sec. Ribbed wall con- 
figuration: 

A = 0.7 in. 
B = 0.125 in. 
TH = 0,06 


N = 10 
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TABLE V 


COMPARISON OF SMOOTH WALL AND RIBBED WALL 
DIFFUSOR EFFICIENCIES* 


SOOM | ----2a th eee 
220 81.21 78.29 0.9640 
rr 69.79 Jo  -.«--- 
WES 65.30 61.75 0.9456 
205 58.58 59.63 1.0179 
30° 53.51 58 44 1.1015 
32° 47 84 56.95 1.1906 | 
34° 43.51 34.55 0.7941 
36° 39.30 29.82 0.7588 
38° 34.55 > if | 45a | 
49° 31.66 27.69 0.8746 


“All data recorded at a measured flow rate of 
approximately 1.78 1b,/sec. Ribbed wall con- 
figuration: 


A = 0.7 in. 

Be = O4425° 1h: 
TH = 0.06 

N = 10 
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4.1.3 Effect of varying the number of grooves upon ribbed wall 


diffusor performance. Table VI is the result of the third and final 
phase of this investigation. Unfortunately, time did not permit this 
phase of the testing to proceed beyond N = 10 grooves. Qne interesting 
result of these tests is illustrated in Table VI. At [is = 30° and 

[3 = 32°, it was noted that the efficiency of the ribbed wall dif- 
fusor decreased by approximately 20 per cent as the number of grooves 
was increased from 7 to 8. This phenomenon did not occur at /3 = 30°. 
Migay, Ref. 1, reports a similar drop in efficiency for a three di- 
mensional conical diffusor operating at an opening angle of 31 degrees. 
His explanation for this decrease in efficiency is that it is caused 

by "'acoustic irradiation of inner-fin hollows.'' No attempt was made in 
this study to either explain or further investigate the phenomenon. 

4.2 Recommendations. Much experimental work remains to be done 
in the investigation of ribbed walled diffusors. It is hoped that 
continued studies in this field will make available data upon which a 
theoretical analysis of ribbed wall diffusor performance can be based. 

It is felt that the next step in the investigation of ribbed 
wall diffusors should be an attempt to obtain data reflecting the 
effect on performance as other parameters of the rib configuration 
are varied. 

In conjunction with further experimental studies into the effect 
of ribs on the performance of diffusors, it is felt that studies into 
the mechanism of flow past grooves or ribs should be thoroughly in- 
vestigated. Fox, Ref. 10, has done much work in this field. It is 
felt that the results of his study would be of great value in attempt- 


ing to understand the mechanism of the ribbed wall diffusor operation. 
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TABLE VI 


EFFECT OF VARYING THE NUMBER OF GROOVES UPON THE PERFORMANCE 
OF A RIBBED WALL DIFFUSOR,.* 


o—— h 2007, S =80ri25, TH =-OvG6 





£ a Bieg) 097, 66 = 05185, TH = o806 


or a ee 





ae A = 0.7, B = 0.125, TH = 0.06 






*All runs for this test made with measured flow rate of 
approximately 1.78 1b /sec. 
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The possible use of ribs, or grooves, as a means of affecting 
boundary layer control ina variety of external as well as internal 
flow situations seems to indicate that further research into the 
behavior of ribs is warranted. It is hoped that the equipment, pro- 
cedures, and experience assembled in the course of this study will 


serve some useful purpose in any such subsequent investigations. 
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APPENDIX A 
COMPUTER PROGRAM 


General statements concerning program "Diffusor.'’ All of the 
calculations discussed in Chapter III were programmed for the IBM 
360-67 Digital Computer in operation at the Computer Facility, Naval 
Postgraduate School, Monterey. The program was written in Fortran IV 
language. All of the calculations as discussed in Chapter III are 
utilized in the program named "'Diffusor.'' Appendix B is a listing of 
the program. Included in the program listing are ''comment'' cards to 
help explain each step of the calculation. A basic understanding of 
the Fortran IV language should be all that is required to follow the 
program operation. In this Appendix the emphasis is directed toward 


the input and output phase of the program operation. 


Program input. The special listing of symbols entitled ''Program 
Input Variables!’ included with this Appendix defines the Fortran 
variables names used as input for program "'Diffusor.'' Figure 13 is a 
sample data page illustrating the order and variable names as recorded 
during the tests. Figure 14 shows the field specification for each of 
the input variables. 

The first data card of each program run is a 72 space "'tiollerith'"' 
Statement. The first 72 spaces of this card may be utilized to include 
supplementary or amplifying statements to be printed as part of the 
program output. Any information punched on this card will be repro- 


duced directly beneath the title line on the output page. 
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MM 
IRUN 
MO 
IDAY 
PATM 
PINOZ 
DLTP 
TN 
TEMP 
PI 
QI 
TTI 
BETA 
TER 
AA 
BB 
TH 


NN 


PROGRAM INPUT VARIABLES 


Number of sets of data to be read 
Number of run 

Month (numeric) 

Day of month (numeric) 

Atmospheric pressure 

Static pressure ahead of nozzle ( P,) 
Pressure drop across nozzle (AP ) 
Total temperature at nozzle 
Temperature (°C) for PATM correction 
Static pressure at inlet Py 

Dynamic head at inlet (P,-P1) 

Total temperature at inlet (Tey) 
Total diffusor opening angle 

Outside ambient air temperature (°F) 
Depth of grooves 

Width of grooves 

Spacing between grooves 


Number of grooves 
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SAMPLE INPUT/OUTPUT D&TA FOR PROGRAM DIFFUSOR 


72 Space Hollerith Statement Card 














ERD Se 







40 
es Be 









FIGURE 14 


INPUT DATA FIELD AND FORMAT SPECIFICATIONS FOR PROGRAM DIFFUSOR 
(Top line indicates card spacing) 
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The second data card of each computer run will indicate the 
number of data sets which are to be read. The name given to this 
"Fixed point'’ variable is MM. 

The remaining data cards represent the input to the program. 
Each data point requires three cards. As illustrated in Fig. 14, the 
first of each series of three input cards contains the following: 
IRUN, MO, IDAY, PATM, PINOZ, DLTP, TN, and TEMP. The second of each 
series of three cards includes: PI, QI, TTI, BETA, and TPR. The last 
in the set of three cards used for each data point contains the fol- 
lowing information related to the ribs of the diffusor: AA, BB, TH, 
and NN. The symbols used to identify the input variables are defined 
in this Appendix. For smooth wall runs, this card was left blank 
resulting in zeros for all of the rib parameters. 

Special attention is directed to the fact that the ''I'' format 
variables are right justified within their specified fields whereas 
the ''F format variables can be punched anywhere within the specified 
field provided the decimal is punched with the number. The field 
specifications and the format for each of the input variables is 
included in Fig. 14. 

Program output. Figure 15 is a sample of the output from program 
"Diffusor.'' The same data used to illustrate the input in Figs. 13 
and 14 are printed as output in Fig. 15. 

The second line of output is the Hollerith statement read into 
the program as the first data card. 

The grouping labeled DATA CHECK is a print out of the raw data 
as read into the program. It was used to check the output for 


possible errors in card punching. 
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PTGURE 15 


SAMPLE OUTPUT PROGRAM DIFFUSOR 


The input sections under the headings of ''F low Rate'' and 


''0;ffusor Performance'' are the converted values of the input data. 
The output sections under these two headings are those values 


calculated in accordance with the discussion in Chapter III. 
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